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0\ \ 1. Introduction 

On 

Several studies have addressed the topic of spiral interactions in pairs or in a 
small group (Barnes 1998). However there is no specific study, to our knowledge, 
of interacting triplets of spirals. Here we present a preliminary report of an 
ongoing research on the dynamics of triplets that includes galaxies of different 
CO \ morphology. In particular we follow a triplet of equal-mass spirals up to a merger 

remnant that resembles an elliptical-like galaxy has been formed. 

It is known that the HI component of spirals is a good tracer of interactions, 
t/^ ■ Indeed, several triplets (e.g. in M81 group) that do not appear to be perturbed 

in optical images do so in HI maps. Although a self-consistent treatment of gas 
requires the inclusion of e.g. pressure forces, we use here a test-particle approach 
to study it during a triple interaction of spirals; this can yield information about 
| the large-scale kinematical features of gaseous tails. 
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2. The Numerical Experiment 

O 1 We use Galaxy-like spirals with total mass M « 5.8 x lO n M , i?haio ~ 135 

£j ! kpc, and disc mass M d w 4.4 x 10 10 M (Kuijken & Dubinski 1995, Model 

B). For computational economy, the number of particles used for each spiral 
were: iVb u i ge = 1000, iVdi sc — 5000 and A^alo — 12000. In isolation the spiral 
was relatively stable for a few orbital periods. Discs' orientations were obtained 
from three random sets of Euler angles (some animated gifs of the simulations 

i ' are at |http : //www . iaa . es/ae/triplet s . html ) . 
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The initial positions and bulk velocities of galaxies were obtained from one 
cold-collapse simulation of spherical galaxy triplets, in which a relatively rapid 
triple merger was found (Aceves 1999). We started the simulation at to w 2 
Gyr from turn- around, when two of the galaxies Gl & G2 had their halos just 
overlapping. These galaxies were approaching at a relative velocity of ~ 140 
km/s; the third galaxy's (G3) halo had not yet 'touch' the binary's common 
halo and was infalling at ~ 40 km/s. In the following, times quoted will be the 
elapsed time since this to. The simulation lasted for another « 12 Gyr. 

The 'gas' particles were distributed uniformly in the plane of the discs 
extending twice the stellar component. Forces on them were calculated using 
the mass interior to their position, with respect to each galaxy centre. The 
distribution of mass in galaxies was assumed to remain the same throughout the 
simulation and that it had a radial dependence. We consider this a somewhat 
better approximation than assuming that all the mass of a galaxy resides at its 
centre; no important extra computational effort was required in this approach. 
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Figure 1. (/e/t) Xy-projection of disc particles of galaxies Gl & G2 
at « 4 Gyr from to- (right) Gas particles of galaxies, as well as the 
bulges, are shown. Numbers are intended to follow to each galaxy. 

3. Results 

In Fig. m we show a projection of galaxies Gl & G2 discs interacting at ^ 4 Gyr; 
galaxy G3 has not yet come close enough for its halo to be perturbed. The discs 
show large damages and a long stellar tidal tail has developed in G2 (~ 150 kpc, 
in projection). The gaseous tails generated reach ~ 200 kpc in projection. Some 
stars and gas particles are returning to the 'binary' galaxy due to conservation 
of angular momentum and because they do not have enough kinetic energy to 
leave the local potential well of the binary. At « 10 Gyr galaxies Gl & G2 have 
already merged with G3; see Fig. [| Galaxy G3 shows a very large projected 
stellar tidal tail of ~ 350 kpc, which is about the same size as its gaseous one. 
At this time the tidal tail of G2 has almost disappear, although important traces 
of it exist as a more diffuse component. 

The configuration at the end of the simulation is shown in Fig. [3[ Enormous 
gaseous tails ~ 1 Mpc are present and the stellar tidal tail of G3 seen in Fig. [2] 
is still present. Since the gas thermodynamics has not been considered, it is 
probable that some 'HI gas' has been ionized by shocks developed during the 
interaction and, hence, this extension can be an upper limit; also, the initial size 
of the gaseous discs and the spiral matter profile determine this extent. After 
removing all unbound particles the shape of the merger was found to be prolate. 
No disc or bulge particles were formally found to be unbound, but galaxies Gl 
& G2 had lost « 3% of their halo mass and G3 « 10%. The merger remnant 
shows a luminous profile that resembles that of an elliptical, although with a 
rather isothermal core, its phase-space has retained more information about its 
spiral past; see Fig. |j. The physical size of the stellar tails reach ~ 500 kpc! 

Figure [5] shows velocity profiles for the gas, which would mimic the result 
of placing a narrow slit along the X-axis, both at « 4 Gyr (top) and 12 Gyr 
(bottom). Velocities along the l.o.s. of ~ 100 km/s in the outer parts with a 
width of ^ (50,100) km/s are obtained, respectively, from Gl & G2 at w 4 
Gyr. At the end of the simulation Gl shows a broad dispersion in the velocities, 
while galaxies G2 h G3 show a more well defined width of ~ 50 km/s in the 
XV^-plane. Tidal tails in G3 developed very late in the merging process. 
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Figure 2. Similar as in Fig. [I], but at « 10 Gyr. Note the change of 
scale. The 'gas' tail is just separating from the stellar one in G3, as 
seen in some observations. The bulge of G3 is shown on the right. 




Figure 3. Final configuration. Note again the change in scale, and 
how some stars are very far away from the central region. 
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Figure 4. (left) Three orthogonal surface number density profiles of 
luminous matter at t w 12 Gyr. (right) Phase-space portraits of discs. 
Note that signatures of the encounter in the three galaxies persist. 
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Figure 5. Velocity profiles of gas along the X-axis at « 4 Gyr (top), 
and at w 12 Gyr (bottom), for each spiral galaxy. 



4. Final Comments 

Some particular results suggested by the present simulation, with its specific 
cold-collapse initial conditions, tend to show that: 

i. Encounters of spirals can effectively destroy stellar discs in <^ bubble • Large 
projected stellar tidal tails ~ 200 kpc can develop during a binary forma- 
tion, and up to ~ 350 kpc in the 3rd. when this joins in. However tidal tails 
are not eternal, they can fall-back to the remnant of their 'parent' galaxies 
in <s ^Hubble? although important traces of them remain for ^> ^Hubble- 

ii. Diffuse light, due to stars stripped from the discs, is to be expected in 
deep imaging of spiral interactions and mergers. However if this stripping 
occurred more than a stellar evolution time-scale ago one may not detected 
it, but stellar remnants may be in the field. 

Hi. Gaseous tidal tails may extend to ^ 1 Mpc in radius, similar to the extent 
of the common dark halo formed, in a triple interaction in ~ bubble • This 
result suggests that gaseous halos may trace the size of a common dark 
halo developed during galaxy interactions. 

iv. Although the stellar merger remnant shows an elliptical-like profile, it is 
more physically justifiable to look for merger signatures in velocity-space; 
albeit the observational difficulties. 
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